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SPECTRAL VARIATIONS IN EARLY-TYPE GALAXIES AS A FUNCTION OF MASS 



o 
o 
o 

< 
(N 



> 

m 

o 
o 
o 

Oh; 

6 



Kristi Dendy Concannon and James A. Rose 

Department of Physics and Astronomy, CB #3255, University of North Carolina, Chapel Hill, NC 27599 

AND 

Nelson Caldwell 

F. L. Whipple Observatory, Smithsonian Institution, P.O. Box 97, Amado, AZ 85645 
Draft version February 1, 2008 

ABSTRACT 

We report on the strengths of three spectral indicators — Mg2, H/3, and Hn/Fe — in the integrated 
light of a sample of 100 field and cluster E/SO galaxies. The measured indices are sensitive to age and/or 
and metallicity variations within the galaxy sample. Using linear regression analysis for data with 
non-uniform errors, we determine the intrinsic scatter present among the spectral indices of our galaxy 
sample as a function of internal velocity dispersion. Our analysis indicates that there is significantly 
more intrinsic scatter in the two Balmer line indices than in the Mg2 index, indicating that the Balmer 
indices provide more dynamic range in determining the age of a stellar population than does the Mg2 
index. Furthermore, the scatter is much larger for the low velocity dispersion galaxies, indicating that 
star formation has occurred more recently in the lower mass galaxies. 

Subject headings: galaxies: abundances — galaxies:stellar content 



1. INTRODUCTION 

An enormous effort has been made in recent years to 
disentangle the degenerate effects of age and metallicity 
on the integrated light of a galaxy. This age-metallicity 
degeneracy (Worthey 1994, 1996) plagues both the pho- 
tometric colors and the integrated spectrum of a galaxy, 
making it extremely difficult to separate the effects of the 
two. Because broadband colors cannot break the degener- 
acy (Bica et al. 1990; Chariot & Silk 1994), the best way to 
determine a unique solution to the problem is to use a com- 
bination of population synthesis models and spectral line 
indices. Absolute ages and metallicities determined using 
this method are naturally subject to systematic errors, yet 
relative values are still of great interest and are generally 
more reliable (e.g., Trager et al. 2000a). Although no sin- 
gle spectral index can completely separate the effects of 
age and metallicity, it is possible to determine the relative 
range empirically by plotting a single spectral index with 
a predominant age or metallicity sensitivity versus a mass 
indicative parameter such as the absolute magnitude or 
internal stellar velocity dispersion. The spread in index 
values which is observed as a function of mass is then di- 
rectly related to the range in age and/or metallicity in the 
galaxy sample. 

The most widely used spectral index system is the set 
of indices originally defined in Burstein et al. (1984) and 
further refined in Worthey et al. (1994) and Trager et al. 
(1998), now commonly known as the Lick system. In this 
system, the H/3 index is used as the primary age-sensitive 
spectral indicator, whereas the Mg and Fe indices are used 
as the primary metallicity indicators. Using the Lick in- 
dices, many groups have been able to place constraints 
on the ages and metallicities of various galaxy samples. 
For example, Trager et al. (2000a) find that for a sample 
of early-type galaxies in low density environments there 



is a large spread in H/3 values (i.e., age), but very little 
variation in metallicity. For galaxies in the Fornax clus- 
ter, though, Kuntschner (2000) finds the opposite effect, 
in that a large spread in metallicity is present with little 
variation in age. Whether such contrasting results are a 
product of the different environments in the two samples 
is an interesting and still unanswered question. It is clear, 
however, that a comparison of the range in age and metal- 
licity of galaxies in different environments will lead to a 
better understanding of the parameters which govern the 
evolution of galaxies. 

To constrain galaxy evolution scenarios, it is imperative 
that we first derive reliable values for the ages and metal- 
licities of galaxies. In this paper we present first results of 
a sample of 100 early-type galaxies covering a large base- 
line in mass and a large range in wavelength in order to use 
both the Lick/IDS age/metallicity indicators as well as ad- 
ditional indices centered around the higher order Balmer 
lines. Although the data will eventually be used to de- 
termine absolute ages and metallicities, and thereby place 
constraints on the theories governing galaxy formation and 
evolution, these goals are beyond the scope of this paper. 
In this letter, we demonstrate the importance of the lower 
mass galaxies and higher order Balmer lines for obtaining 
a complete view of the ages of early- type galaxies. The 
paper is organized such that the data and reduction pro- 
cedures are summarized in §2, the spectral indices and the 
index diagrams are presented in §3, and the analysis of 
these diagrams is given in §4 and discussed in §5. 

2. OBSERVATIONS 

Our galaxy sample consists of 100 early-type galaxies, 
70 of which are field galaxies and 30 of which are in the 
Virgo cluster. The field galaxies are drawn from the CfA 
redshift survey (Huchra et al. 1983); the cluster galaxies 
are from the Virgo Cluster Catalog (VCC) of Binggeli et 
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al. (1985). The observed galaxies cover a range in absolute 
magnitude of — 16 < Mb < —22 and have been previously 
classified as either E or SO. 

Long slit spectra of the galaxies were obtained at the 
F. L. Whipple Observatory with the 60" Tillinghast tele- 
scope, the FAST spectrograph and a Loral 512 x 2688 
pixel CCD (Fabricant et al. 1998) during seven different 
observing runs between 1997 September and 1999 Octo- 
ber in seeing conditions of 1" or better. A 3" slit and 
the 600 line/mm grating were used to give a dispersion of 
0.75A/pixcl and 3 A FWHM spectral resolution and a spec- 
tral coverage of 3500A - 5500A. Each exposure was 1800 
seconds with 2-4 exposiires taken for each galaxy. The slit 
was oriented at the parallactic angle so as to preserve the 
blue light, resulting in a S/N per pixel between 70 and 100 
at 4000A in the final coadded spectrum. 

The spectra were reduced in IRAF with standard re- 
duction procedures including bias and dark subtraction as 
well as flat fielding. The spectra were wavelength cali- 
brated using an internal HeNeAr lamp and were flux cali- 
brated using standard stars from the Massey et al. (1988) 
spectrophotometric catalog. Multiple galaxy exposures 
were added together and all galaxies were de-redshiftcd 
to the same rest wavelength. Velocity dispersions were de- 
termined using the Fourier cross-correlation method (e.g., 
Tonry & Davis 1979) in the task fxcor and various stellar 
templates which were observed during twilight. All spec- 
tra were then smoothed to the same effective dispersion of 
230 km/s to insure that each object experienced the same 
intrinsic broadening. 

3. SPECTRAL INDICES 

Two types of spectral indices were measured from the 
calibrated and smoothed spectra, the pseudo-equivalent 
width indices of the Lick group (Burstein et al. 1984; 
Worthey et al. 1994) and the line ratios of Rose (1994; 
CaldwcU & Rose 1998). The Lick/IDS family of indices 
measiircs the pseudo-equivalent width of absorption line 
features such as H/3 and the Mg bands which lie redward 
of 4000A. The Lick equivalent width indices are measured 
by centering a bandpass of typically 30-40A on the absorp- 
tion feature and using a red and blue sideband of similar 
width as the continuum reference. In addition to the Lick 
system, we have also chosen to work at bluer wavelengths 
where young stars contribute a larger fraction of the light. 
The indices used in this region are the Rose spectral line 
ratios, which arc defined as the ratio of the residual in- 
tensity in two neighboring spectral lines. For our galaxy 
sample, we have measured a total of 36 spectral indices, 22 
Lick indices and 14 Rose indices, and report in this letter 
on three of the key age/metallicity diagnostics. The spec- 
tral indices analyzed here are: Mg2, a Lick index which 
measures the strength of the Mg H -|- Mg b molecular 
bands, H/3, another Lick index which measures the equiv- 
alent width of the H/? feature at 486lA; and Hn/Fc, a Rose 
spectral line index which is an average of three hydrogen 
to iron line ratios^ and is defined such that a lower index 
value corresponds to a stronger Balmer line strength. 

For each index, the error for an individual galaxy is 
determined from the rms scatter in the repeat measure- 
ments. Since this method is unreliable for those galaxies 

^Hn/Fe =< H5/Fe4045 + H'y/Fe4325 + H8/Fe3859 > 



which have only two observations, we have defined our er- 
rors for a particular index and galaxy utilizing the errors 
in all indices and galaxies in the following maimer. The 
error, Aj^ , in index j for galaxy i (where there are a total 
of m indices and n galaxies) is given by: 

En 
_ i=i Sjj 2^j=i Sij 

where ^"-i £ij is the rms error in the j*'' index averaged 
over all galaxies, Y^^i ^ij is the error in the i*^ galaxy av- 
eraged over all indices of the same type (i.e, all Lick indices 
or all Rose indices), and "Y^^'j^i ^ij is the average error in 
all indices of the same type for all galaxies. The aver- 
age fractional uncertainty in each index is 1% for Hn/Fe 
and Mg2, 3% for H/3 and 2% for the velocity dispersion 
measurements. 

The three indices reported in this paper are shown in 
Figure 1 as a function of log a (i.e., as a function of mass). 
The line is the best fit to the high-mass data extended 
to the low-mass end; an explanation of this fit is given 
in §4. Comparing these diagrams, two effects are immedi- 
ately evident. First, the scatter among galaxies appears to 
increase between the high-mass end and the low-mass end 
of the diagrams. Second, this increase in scatter is signifi- 
cantly more dramatic in both the H/3 and Hn/Fc diagrams 
than it is in the Mg2 diagram. Some of the observed scat- 
ter of the galaxies about the linear fits is obviously due 
to measurement errors associated with the data, but vari- 
ations in the age and mctallicity of the galaxies play a 
contributing role as well, and in practice, it is difficult to 
separate this intrinsic scatter from the observational er- 
rors. In the following section, we describe a statistical 
analysis of the scatter in the index diagrams and discuss 
their implications. 

4. RESULTS 

In each diagram in Figure 1, a linear relation is readily 
identifiable for the galaxies with log a > 2.0. This prop- 
erty is related to the well-known color-magnitude relation 
of elliptical galaxies (Visvanathan & Sandage 1977) and is 
thought to be driven by metaflicity (Faber 1973, 1977). 
The relation is naturally explained by the supernova- 
driven wind model (Larson 1974; Arimoto & Yoshii 1987) 
in which the more massive galaxies are better able to re- 
tain their supernova ejecta than the smaller galaxies and 
thus become more metal-rich and redder. Because the 
high-mass galaxies represent a fiducial, normal sample of 
galaxies, we used only these points to determine the best 
fit to the data and extended the fit to the lower mass galax- 
ies. The three outlying points above log a = 2.0 were not 
used in the regression analysis, since each of these galax- 
ies is known to be in an interacting system with star for- 
mation, possibly triggered by this interaction. The least- 
squares fitting is done using the BCES (bivariate corre- 
lated errors and intrinsic scatter) method of Akritas & 
Bershady (1996) that fits data for which there is some in- 
trinsic scatter present and for which the error bars are 
heteroscedastic (non-uniform) in both the x and y coordi- 
nates. Because we are interested in measuring the intrinsic 
scatter as a function of mass (i.e., internal velocity disper- 
sion), the chosen fit minimizes the residuals in the indices 
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(?/- variables), which have, on average, larger errors than 
the velocity dispersion (x-variable). To verify the fits, the 
analysis was repeated for 1000 bootstrap re-samplings of 
the data. We find that the BCES and bootstrap fits agree 
within their associated errors. 

To separate the intrinsic scatter, rj, from the scatter due 
to measurement errors, we used an estimate given by Fioc 
& Rocca-Volmerange (1999): 

= Y^lbn Bx, Af {B^e, + eyd] (2) 
where B and A arc the slope and intercept respectively 
of the BCES regression analysis, Xi is the observed log cr. 
■tji is the observed index value and ^xi, iyi are the asso- 
ciated errors in the Xi and yt observations. The galaxies 
were divided into two bins, log cr > 2.0 and log a < 2.0, the 
scatter was calculated for each subsample, and the analysis 
was verified with Monte-Carlo simulations. The measured 
intrinsic scatter for both the high- and low-mass subsam- 
ples is shown as the vertical dotted lines in the lower right 
corners of Figure 1 . 

The results of the statistical analysis of the three spec- 
tral indices, H/3, Mg2, and Hn/Fe, are shown in Table 1. 
Here, the intrinsic scatter in each subsample is denoted by 
ry; the error in 77, estimated from bootstrap re-samplings, is 
given by (5^. The reader should be aware that because the 
two Lick indices and the Rose index are defined differently 
(as equivalent widths and a line ratio, respectively) it is 
not meaningful to compare the absolute difference in the 
intrinsic scatter between two indices, i.e., it is only mean- 
ingful to compare relative differences in the same index. 
The quantitative analysis clearly verifies our previous pre- 
dictions in that a non-zero intrinsic scatter exists in each of 
the index diagrams and that this scatter is greater for the 
low-mass galaxies than for the high-mass galaxies. In fact, 
the intrinsic scatter increases by a factor of 1.5 from the 
low-mass to the high mass galaxies in the Mg2 diagram, 
by a factor of 2.6 in the H/3 diagram, and by a factor of 5.2 
in the Hn/Fc diagram (R value in Tabic 1). In addition 
to the three indices reported here, we have also estimated 
the intrinsic scatter in another Lick index, Mg b, which has 
become a popular age/met allicity indicator (e.g, Trager et 
al. 2000a). For this index, we find that the intrinsic scat- 
ter increases by a factor of 2.2 from the high- mass to the 
low-mass galaxies. It is worth noting that this value is 
between the estimates for Mg2 and H/3. 

Prom this analysis, it is evident that more intrinsic scat- 
ter exists in the Balmer indices than in the Mg2 index. If 
the factor driving the scatter were metallicity instead of 
age, we would actually expect to see a larger scatter in the 
Mg2 data than in the Balmer indices, since Mg2 is more 
sensitive to metallicity than to age in comparison with H/3 
and other Balmer lines (e.g., Buzzoni et al. 1994). How- 
ever, the opposite situation is observed. Because we see 
more variations in the Balmer indices than the Mg2 index, 
especially at low mass, the index diagrams imply that the 



lower-mass galaxies exhibit a larger range in ages than the 
high-mass galaxies. 

5. DISCUSSION 

The results of this analysis bring forth two key issues 
that must be addressed. First, since the two Balmer line 
indices exhibit more intrinsic scatter than the Mg2 index, 
this qualifies them as more robust age indicators. The 
relatively little intrinsic scatter observed in the Mg2-log 
cr diagram indicates that this index is not well suited for 
constraining the age of a stellar population. The fact that 
significant age variations may be hidden in the Mg2-log 
(T relation has been previously noted by both Worthey et 
al. (1996) and J0rgensen (1999) and has been studied in 
detail by Trager et al. (2000b). This finding may have a 
significant impact on the previous results of authors who 
have used the Mg2 index as their principal age diagnos- 
tic (e.g., Bernardi et al. 1998), and their conclusions may 
need to be re-evaluated using a more age-sensitive index 
such as H/? or Hn/Fe. 

The second key result is the amount of scatter observed 
for the low-mass galaxies in the index diagrams. The dra- 
matic increase in scatter, particularly at the low-mass end 
of the Balmer line diagrams, indicates that the smaller 
galaxies have experienced a more varied star formation his- 
tory and that these galaxies have a larger range in age than 
their brighter counterparts. This finding places significant 
constraints on a hierarchical formation scenario such as 
the CDM model (Baugh, Cole & Frenk 1996; Kauffmann 
1996), since in this case one would naively expect to see 
that low-mass galaxies are older than the larger galaxies 
(Kauffmann & Chariot 1998). However, because the ma- 
jority of the low-mass galaxies in our sample are in the 
Virgo cluster, it is difficult to determine whether the ob- 
served variation in age is simply a product of the differ- 
ent environments. Even if this were the case, though, in 
a hierarchical scenario we would expect that galaxies in 
clusters would be older than field galaxies, since merging 
will happen more rapidly in clusters. Thus, regardless of 
the environmental influence, our results appear to place 
significant constraints on the hierarchical picture. 

In a forthcoming paper, we will report on use of the 
higher order Balmer index and population synthesis mod- 
els to determine the spread of age and metallicity of our 
sample, and thereby place quantitative constraints on for- 
mation and evolution scenarios. 

We would like to thank to William J. Thompson for 
much valuable advice and discussion concerning the sta- 
tistical analysis in this paper. Thanks are also due to 
Michael Akritas and the SCCA^ for their advice. This re- 
search was supported in part by a North Carolina NASA 
Space Grant Consortium Graduate Fellowship and by NSF 
grant AST-9900720 to the University of North Carofina. 
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Fig. 1. — Three spectral indices (a) Mg2, (b) H/3, (c) Hii/Fe plotted versus log a. The solid circles arc field galaxies; the open circles are 
Virgo cluster galaxies. The regression line is the calculated fit with the dashed portion indicating the region of the line which is extrapolated. 
The vertical dotted lines represent the measured amount of intrinsic scatter in each subsample; the high-mass scatter is on the right, the 
low-mass scatter is on the left. The index values in (a) are in magnitudes; the values in (b) are equivalent widths in Angstroms; the index in 
(c) is a ratio. 
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Table 1 



Measured Intrinsic Scatter in Index Diagrams 



Index 


Subsample'^ 


V 


Sn 






Mg2 


low mass 


0.050 


0.005 


1.52 


0.21 




high mass 


0.033 


0.003 






H/3 


low mass 


0.578 


0.084 


2.60 


0.48 




high mass 


0.222 


0.025 






Hn/Fe 


low mass 


0.105 


0.015 


5.25 


1.01 


high mass 


0.020 


0.003 







''Low mass: log a < 2.0; high mass: log a > 2.0 

Ratio of intrinsic scatter in low mass galaxies to that 
in high mass galaxies 

'^Error in the intrinsic scatter ratio 



